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Abstract 

We give a constructive proof for the superbosonization formula for 
invariant random matrix ensembles, which is the supersymmetry ana- 
log of the theory of Wishart matrices. Formulas are given for unitary, 
orthogonal and symplectic symmetry, but worked out explicitly only 
for the orthogonal case. The method promises to become a powerful 
tool for investigating the universality of spectral correlation functions 
for a broad class of random matrix ensembles of non-Gaussian type. 

Pacs: 0250.-r, 0540.-a 



1 Introduction 

I will report on superbosonization, a technique in random matrix theory 
(RMT), the appellation given by Efetov et al. [T] which we want to make a 
rigorous tool. Recent mathematical and theoretical work from SFB/TR12 is 
to be published [21 13] . The method is the supersymmetry analog of the theory 
of Wishart matrices. Other people have done related work connected with 
universality properties of invariant ensembles [U El El E]. But nowhere so 
far the precise integration domain for the supermatrices under consideration 
have been specified. Here I want to present a constructive proof of the 
superbosonization formula on a level, which can be understood by a physicist. 



* Presented at the Conference on Random Matrix Theory: From Fundamental Physics 
to Applications, Krakow, Poland, May 2-6, 2007 
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One may generate correlations of random Green functions from some 
Gaussian integrals over a set of A'"-dimensional vectors building a, N x p 
dimensional matrix which after averaging are no longer Gaussian, but 
the integrand is a function of the matrix Q = of invariants. The 

transition to this new set of variables is called bosonization and it makes e.g. 
large- evaluations simpler using saddle-point techniques. This kind of order 
parameter was previously introduced in the theory of spin glasses and the 
theory of Anderson localization. These theories involve the socalled rephca 
trick, where at the end of calculation, after averaging over the ensemble, one 
has to take the limit of number of replicas to zero, a procedure that is very 
difficult to make mathematically rigorous. 

A more rigorous method is to introduce supervectors ^ (if possible) with 
commuting and anticommuting components. Then no replica limit is needed, 
in some sense the dimension of Grassmannian (anticommuting) vectors can 
be considered as negative. The change of variables from a set of supervectors 
\I' to a supermatrix Q = j^'if^'if is called superbosonization. It is useful for 
invariant ensembles and non-Gaussian ensembles and avoids the sometimes 
cumbersome Hubbard-Stratonovic transformation. 

2 Superbosonization formula 

Let us assume, that we want to integrate a function F of invariants (with 
respect to a symmetry group) 



Here C are anticommuting variables building a rectangular N x p matrix (p 
even), $ arc commuting variables building a, N x q matrix (and we need 
lateron N ^ q). Originally one starts with an invariant function of C, and 
$ and ends up with a function of invariants, which is not unique for Grass- 
mannians, nevertheless the result of integration is unique. DC, is the Berezin 
integration form and is the flat measure of matrix elements. 

We want to transform to integrations over supermatrices Q = ( ^ ) with 
bosonic (commuting) entries A, B and fermionic (anticommuting) entries cr, 
(7^ . In this way we reduce considerably the number of bosonic and fermionic 
integration variables especially in the case where N goes very large. It is not 
surprising that B runs over positive Hermitian matrices and a and cr^ run 
over Grassmannians. But the amazing thing is that the entry in the Fermi- 
Fermi sector, which is originally nilpotcnt, is replaced by a matrix A which 
runs over a manifold of unitary matrices. There is a compact way of writing 
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the integral ([T]) as superintegral containing some power M of the socalled 
superdeterminant S det Q of the supermatrix Q: 



•■pq 



M DA DB D(a,a^ 



det 5 



H M 



F 



A at 
a B 



_det{A - a^^a) 

= M j DQ{^ det QY^ F{Q) (2) 

The appearing measures are the flat measures and M is a, normalization 
constant. This is the superbosonization formula and I want to specify it for 
different symmetry groups. 

For the real orthogonal group On the dagger ^ just means the trans- 
posed A is skew-symmetric {A = —A^) and unitary, and B is real symmet- 
ric {B = B^) and positive. The power M is given by M = {N + p — q — l)/2. 
For the unitary symplectic group USp^ (-^ even) the dagger ^ means the 
dual ^, with e.g. A^ = ZA'^Z'^ where Z is the symplectic unit which is 
quasidiagonal with quaternion elements ( i*! o ) the diagonal. A is anti- 
selfdual [A = —A^) and unitary and B is Hermitian, selfdual {B = B^) and 
positive. The power M is given by M = (A^ + p — q + l)/2. In the case of 
unitary group Uiy the dagger means just the usual adjoint = transposed and 
complex conjugate. The Grassmannian cjt can be any independent set a of 
Grassmannians. Here M = N + p — q. 

Let us say a bit on the manifolds. A = —A^ implies that (ZA)^ = 
{ZA), i.e. for orthogonal symmetry ZA belongs to the set of selfdual unitary 
matrices, which is called in physics the circular symplectic ensemble CSE 
(not only meant as an invariant measure but also as a manifold). Similarly 
A = —A^ implies that ZA is symmetric. Thus for symplectic symmetry 
ZA belongs to the circular orthogonal ensemble COE. Finally for unitary 
symmetry ZA belongs to the circular unitary ensemble CUE. 

On the other hand the matrices B are positive Hermitian and symmetric 
or selfdual in the orthogonal or symplectic case. All those manifolds are 
Riemannian symmetric spaces, which were under special consideration in 
our SFB/TR12, related to the socalled ten-fold way of universalities in RMT 
for fermionic systems [H]. Moreover the supermatrices Q = ( f i)Q belong 
to Riemannian symmetric superspaces introduced by Zirnbauer [9J. 

The measures which we have indicated so far are the flat ones and they 
are related to the invariant measures on the corresponding manifolds by some 
power of determinant or superdeterminant. E.g. the invariant measures on 
the supermanifolds are (up to normalization) 

dj^iQ) = DA DB D{a, a^) (S det Q)^ (3) 
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with R = p- q for (3 = 2, R={p-q - l)/2 for f3 = 1, R = {p - q + l)/2 for 
(3 = 4. The invariant measures on the manifolds of A and B can be read off 
from these expressions by putting q = or p = 0. 

Let us say a httle bit about the normahzation constant 

Af = C{N,q)-D{N-q,p) (4) 

It factorizes in two constants coming from the Bose-Bose sector C {N, q) and 
from the Fermi-Fermi sector D{N — q,p). They are essentially ratios of 
volumes of the corresponding symmetry groups. They are separately only 
defined for N ^ q because N — q is the smallest dimension of symmetry group 
that appears. And the proof shows that this is actually needed. However, 
the product M = C ■ D is even defined for A^ + p — g^O, itis not clear if 
this has physical relevance. 



3 Idea of proof 

Let me now give an idea of the proof of the superbosonization formula. We 
start withp = 0, no Grassmannians, this is just the case of Wishart matrices. 
In this case one has to integrate some function over bosonic vectors 

j D^H{^^^) = j DBD^ 5{B - ^^^)H{B) 

= C{N,q) j DBH{B){detB)^ (5) 

The last equation can simply be found by rescaling $ = with B > 

for N ^ q. Since C{N,q) is independent of the function H{B) we may 
choose H{B) = exp(— Tr B) and find for example with (3=1 using Selberg's 
integral [10] and diagonalization of B with orthogonal matrices: 

B>0 

= ^n^''^'r((iV-g + l+j)/2). (6) 
i=o 

This constant is related to the ratio of volumes of orthogonal groups 
V{ON-q)/V{ON)- Since the manifold i? > is noncompact there appear 
noncompact integrals over eigenvalues of B: bi > 0. 



j DBe-^'^'detB^''-'^-'^/' 



4 



More interesting is the case g = 0, only Grassmannians, the case of 
Grassmann Wishart matrices (: {N — q) x p, p even. In this case the naive 
introduction of a (5-function for A = leads to inconsistencies and we have 
to be more careful. We have to integrate a function G{('^() which we write 
with the help of a shift operator 

DCGieO = / /^C exp f Tr C^C A V(^) 



6A^ 

'det(2^) G{A) 



A=0 



(7) 



A=0 

Here A is an antisymmetric matrix and the ^-integration yields the (A^ — g)-th 
power of the Pfaffian ^det (2^). 

The superbosonization formula to be proven says: 

/DA 
(det A) 2 

where A runs over the unitary antisymmetric matrices. G{A) may be any 
polynomial or analytic function. Since A is unitary there is no problem 
at A = or det A = for this integral. Obviously it is enough to prove 
this formula for any exponential function G{A) = exp(Tr AB/2) where B 
is again antisymmetric. In this case the Pfaffian action ([7]) can simply be 
calculated and the last equation can again be proven by rescaling choosing 
ZB a selfdual unitary matrix. To calculate the constant D{N — q,p) one may 
choose B = Z. Then one finds diagonalizing ZA with symplectic matrices 



DiN-q,p) ViUSp,)p/^{l)\ 



]^(aj - a^yJ^rfafcO;. ^ ' ■ 'e" 

i<j k=l 



2m y/'f^ vr^V^ 



2^~V |ir((iV-g+l+j)/2)' 

This is again related to the ratio of volumes of real orthogonal groups 
V{0]^_g+p)/V{ON_g). Since the manifold ZA G CSE is compact there ap- 
pear integrals over eigenvalues of ZA along the unit circle. Interestingly, 
these compact integrals in the complex plane are related to the noncompact 
integrals with bi > along the real axis. 

Finally let me consider the full supermatrix Q. I only report here on the 
orthogonal case. In the integral 
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we have to go from N x {p, q) variables (C, $) to the (p, q) x (p, q) supermatrix 
Q. We choose an orthogonal matrix O which rotates $ to a quadratic qx q 
matrix \fB 

^ = 0%-%=( 1 . (11) 



;i2) 



B 

Then the Grassmannians are rotated correspondingly 
which has a q x p component Co which transforms to a 



(13) 



and a remaining (A^ — q)xp component (i which can be integrated out using 
(El) and 

C^C = C^C = CiCi + Co^Co . (14) 

Similarly the integration over B follows from ([5]). The result is almost what 
we want 

{N-q+p-l)/2 



Ipg=Af I DA I DB I Da 



detB 
detA 



A + a^^a 



a 



B 



(15) 

Now shift A A — a'^^a, which is possible since the manifold CSE has no 



boundary (compact symmetric space) and we end up with 
Ip, = AfJ DQ (Sdetg)(^-'^+^'-^)/2 ^(g) 



(16) 



with SdetQ = deti?/det(y4 — a^a). This completes the proof for the 
orthogonal case. Similar considerations lead to the corresponding results for 
the symplectic and unitary cases. 
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